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We present a general analysis of the state obtained by subjecting the output from a continuous- 
wave (cw) Gaussian field to non-Gaussian measurements. The generic multimode state of cw Gaus- 
sian fields is characterized by an infinite dimensional covariance matrix involving the noise corre- 
lations of the source. Our theory extracts the information relevant for detection within specific 
temporal output modes from these correlation functions . The formalism is applied to schemes for 
production of non-classical light states from a squeezed beam of light. 
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I. INTRODUCTION 

Quantum computing and communication, and quan- 
tum metrology with atoms and light are topics of high 
current interest offering on the one hand proposals for 
new revolutionary technologies and addressing on the 
other hand fundamental questions about the properties 
of quantum systems. Light is an ideal carrier of infor- 
mation over long distances, and it is an ideal physical 
probe for the properties of matter. Numerous proposals 
have been made for the use of light for quantum infor- 
mation purposes and both theory and experiments have 
advanced to the quest for and application of specially 
designed quantum states of light. The observation that 
coherent states of light suffice for quantum cryptogra- 
phy |]J has not removed the interest in single-photon and 
other non-classical states, since it has been proven that 
the entire family of Gaussian states (to which the coher- 
ent states belong) cannot be distilled by Gaussian oper- 
ations 0, and since, e.g., the implementation of qubits 
for digital quantum computing also necessarily involves 
non-Gaussian states 0- 

Non-Gaussian states can be produced by single photon 
emitters such as atoms or ions in cavities 0, || , or N- V 
centers (for a recent review of single-photon emitters 
see 0), but they may also result from post selecting the 
output of Gaussian light sources after non-Gaussian mea- 
surements have been carried out on part of the field @ . 
Thermal, coherent and squeezed light are all examples of 
Gaussian states, and prompted by recent experimental 
and theoretical work on single photon production from 
a cw squeezed light field @, ES El 13 111 13 , we shall 
in this paper develop a theory for the state of a specific 
output mode, obtained from a cw Gaussian state light 
after an arbitrary measurement has been performed on 
a designated trigger mode. This analysis in fact consists 
in two straightforward steps: 

1. From the general multimode Gaussian state, fully 
characterized by the correlation functions of the light 
beam, we derive the covariance matrix and correspond- 
ing Gaussian Wigner function for the state of the (trig- 



ger+output) two mode system prior to detection. 
2. The two-mode Wigner function has two complex (four 
real) arguments, and the effect of the measurement on 
the trigger mode is described by a (generalized) projec- 
tion, which amounts to a specific operation on the Wigner 
function followed by an integration over the trigger mode 
variables. 

These steps, which can both be carried out analytically 
for several relevant cases, produce a Wigner function for 
the quantum state of the output mode which provides all 
available information about the state. 

In Sec. II, we show explicity how to extract the two- 
mode Gaussian state from the noise correlations of the 
Gaussian field light source. In Sec. Ill, we discuss the 
state reduction and the output field resulting from differ- 
ent kinds of photon counting measurement on the trigger 
mode. In Sec. IV, we specialize to the interesting case 
of single photon production, and we present a concrete 
analysis for the squeezed output from an OPO and we 
summarize the technical implementation of various ele- 
ments and phenomena (beam splitters, losses, inefficient 
detectors, dark counts, ...). Sec. V concludes the paper. 



II. GAUSSIAN STATES - REDUCTION TO 
TWO MODES 

In a generic description of the experiments that we 
wish to analyze, a cw Gaussian light field is split and pos- 
sibly filtered by a sequence of linear optical devices, and 
one particular field mode (in our example: the field at a 
specific instant of time) is registered by a measurement 
device, see Fig. 1. Conditioned on the measurement out- 
put, another component of the beam, and in particular 
the field pertaining to a single designated output mode, 
is emitted in a non-Gaussian state. In the case of a low 
flux from an OPO, the total field contains vacuum and 
a weakly populated two-photon component, and a sin- 
gle detector click, stemming from the two-photon com- 
ponent, suggests that a one photon state is produced. 
At higher flux, the squeezed output is a superposition of 
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Figure 1: Experimental set-up for production of a non- 
Gaussian state from a Gaussian cw source of light. The field is 
split into trigger and output components, and after filtering 
(optional), a non-Gaussian photon-counting detection takes 
place on the trigger mode. Conditioned on the trigger signal, 
the output beam may be in a non-classical state, and by use 
of a temporally modulated local oscillator field, a specific sin- 
gle mode in the output beam can be extracted for homodyne 
detection. 



even photon number states, and a single photon subtrac- 
tion returns a superposition of only odd number states, 
which is a non-classical state with some Schrodinger Cat 
character. This paper addresses how good single pho- 
ton states and how much non-classical character can be 
obtained in a realistic physical set-up. 

Our analysis proceeds by noting that prior to the count 
event, we have a system with two relevant field modes: 
mode 1 : the temporal mode on which the trigger detec- 
tion takes place. We assume an instant click event at 
time t c , but the field may well be frequency filtered prior 
to detection, so that the field registered is actually an 
integral over time of the field that left the original light 
source. 

mode 2: the output mode, populating a quantum state 
that we shall characterize fully by its phase space Wigner 
function. The Wigner function is experimentally avail- 
able by homodyne detection and quantum state tomog- 
raphy |15| , where the homodyne detection can be carried 
out with a time dependent local oscillator field, which 
will extract the desired output mode function. For our 
concrete application we expect the optimum mode func- 
tion to be localized around the time of the click event on 
mode 1, and that it has a temporal shape reflecting the 
band width of the squeezed light source. See Fig. 2. 

Let a(t) and ^(t) denote the annihilation and creation 
operators of the source field. The input to our analysis 
is the noise correlations represented by the correlation 
functions (d(i)d(t')) and (a> (t)a(t')). We do not have 
to restrict the analysis to this case, but we note that 
for a stationary source, these correlation functions only 
depend on the difference between time arguments (t — 
t'). The most general formulation of the trigger+output 
modes and their properties is obtained by designating 
field mode operators, 



Figure 2: Schematic representation of the trigger (1) and out- 
put (2) signals. The OPO produces photons in pairs inside 
a cavity, and when a discrete trigger photon is detected (ver- 
tical lines in the figure), we expect to have a second photon 
in the output beam, temporally separated from the trigger 
event according to a mode function (shaded function in the 
figure) with a width related to the squeezing band width of 
the source (part a). If the squeezed light source is too bright, 
the click events occur more often, and the mode functions of 
the output single photons states start to overlap in an inco- 
herent manner (part b). 



where the functions fi(t) with i = 1,2 denote the mode 
functions of the trigger and output modes respectively, 
taking into account completely arbitrary manipulations 
of the field states by linear filters and beam splitters, 
and detection in arbitrarily shaped mode functions. The 
vacuum fields bi^ vac are needed for unitarity, but their 
effects are fully accounted for by the source terms. 

With explicit expressions for the functions in Q and 
for the source correlations, we can compute all second 
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of the trigger-output modes prior 
to the trigger detection. The cw multi-mode field state 
is Gaussian, and when the state is restricted to only two 
modes (all others are traced out) the Gaussian prop- 
erty remains. This implies that the two-mode Wigner 
function of the field modes 1 and 2 is a Gaussian func- 
tion of four real quadrature variables (2/1,2/2,2/3,2/4) = 
(xi,pi,X2,P2), fully characterized by a covariance ma- 
trix with elements =< yiUj > + < j/jj/j >, expressed 
explicitly in terms of (didj), 



W v (y) = 



1 



TrVdet(V) 



exp(-y T V l y) 



(2) 



= / fi{t')a(t')dt' +{ik 



(1) 



where y is a column vector, and - T denotes the transpose. 
(We assume without loss of generality, that there are no 
mean field components and hence no displacement of the 
state to take care of). 



III. NON-GAUSSIAN STATES 

The Wigner function J5J is so simple, that it consti- 
tutes a natural starting point for analysis of the state of 
the output mode 2, conditioned on various possible mea- 
surements performed on mode 1. Homodyne or hetero- 
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dyne measurements are readily accounted for by elemen- 
tary transformations of the covariance matrix, returning 
a Gaussian state of the output field with a mean ampli- 
tude that depends on the trigger measurement outcome, 
and with a covariance matrix that is independent of the 
actual outcome 0] . The purpose of the present work 
is to go beyond the Gaussian states, and for this purpose 
more general transformations of the Wigner function @ 
must be allowed. Such transformations follow if the trig- 
ger mode is made subject to a measurement that does not 
preserve the Gaussian character, for example a photon 
counting experiment, and we shall analyze three variants 
of photocounting measurements: 



A. Perfect number detector 

A perfect photon number detector will yield an integer 
eigenvalue, n — 0, 1, when applied to the trigger mode. 
Such a detection event is accompanied by a projection of 
the field state on the corresponding trigger mode eigen- 
state with the Wigner function W n (xi,pi), which for the 
remaining output mode implies a change into a state de- 
scribed by the conditioned Wigner function 

W c ,n{x2,P2) =M n J dx\dp\ Wv {x\ , pi , X 2 , P2) W n (xi , pi), 

where the normalization constants Af n are inversely pro- 
portional to the probabilities for the different number 
state outputs. The most interesting case is probably the 
one of a detection of a single photon with 

W n=1 (x u pi) = ^{2x1 + 2p\ - 1) exp(-^ _p2) t (4) 

and we observe that the integral only involves a Gaus- 
sian and a product of quadratic terms and a Gaussian. 
It can hence be carried out analytically, the result being 
itself a product of a Gaussian and a second order poyno- 
mial in the variables X2,P2- Kim et al ^| have studied 
the negativity of the Wigner function at the origin as a 
suitable measure of the non-classicality of the output of 
such a protocol. Their study dealt only with a physical 
two-mode system, but they assumed a sufficiently general 
covariance matrix V that they describe also our general 
cw setting, where V is the result of an integration over 
the cw field correlation functions and arbitrary trigger 
and output mode functions. Their analytical expression 
for the conditioned Wigner function W c ,n=i (^2,^2) and 
in particular for its possible negativity at the origin, may 
be helpful in the search for optimum trigger and out- 
put mode functions in such an experiment. Dakna et al 
[2| have studied the case of a pure single-mode squeezed 
state, which after beam splitting is a special case of the 
Gaussian two-mode states, and they have presented vari- 
ous results including output mode Wigner functions con- 
ditioned on the detection of different number states of 
the trigger mode. With Eas. l|2l3|) . and the appropriate 



covariance matrix, determined from the source correla- 
tion functions, we are able to generalize these results to 
search for Schrodinger cat-like states in cw experiments. 

B. "On/off" detector 

The so-called on/off detector, is able to distinguish be- 
tween a vanishing and a non- vanishing photon number in 
a field mode, but it is unable to resolve photon numbers 
larger than unity. The interesting "on" output of such 
a detector, projects the system on the space spanned 
by all non-vanishing photon number states, equivalent 
to removal of the vacuum component, and the resulting 
Wigner function of the output mode is 

W^ n {x 2l P2) =A/" on J dx 1 dp 1 W v (xi,pi,x 2 ,P2) 

(l-W n=0 {x uPl )), (5) 

where W n= o(x\,pi) = ir^ 1 exp(— x\ — Pi), and the in- 
tegral is readily evaluated for any covariance matrix V. 
The result is readily seen to be the difference between 
two Gaussian functions in the variables X2,P2- This de- 
tector model was recent ly a pplied by Sasaki and Suzuki 
in a theoretical analysis [l^ of the cw OPO output field. 
Their approach did not reduce the system first to the 
two relevant modes, but used instead a mode expansion 
on a complete basis of prolate spheroidal wave functions 
and suitable generating functions for the number distri- 
butions studied also by Zhu and Caves [Isj. They do, 
however, also obtain Wigner functions expressed formally 
as the difference between two Gaussian functions with 
different pre-factors and widths. Given that both the co- 
variance matrix V and the integral JSJ can be evaluated 
analytically, we advocate that our approach provides a 
more straightforward analysis of the "on/off" detector. 

C. Click event 

The most natural and realistic detector model, is the 
one of usual photo detection theory, in which the trigger 
mode is incident on a physical material, from which a 
bound electron is excited into the continuum and sub- 
sequently detected. In this scenario, it is the electron 
that is being detected and made subject to a projection 
into the continuum state space, and the back action on 
the field state is given by the field annihilation opera- 
tor |19| . In a density matrix formulation this amount to 
the transformation p — ► a-\ pa\ /Tr(q| ai p), and using the 
operator correspondence (20j between application of the 
operator (ai = (xi + ip)j\[2 and multiplying and differ- 
entiating the Wigner function with respect to its argu- 
ments x\ , p\ , we find the non-Gaussian two-mode Wigner 
function of the field after the click of the detector. The 
reduced state of the output field is at this point obtained 
by tracing over the trigger mode field, i.e., by integrating 
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the Wigner function with respect to the xi,p± arguments 
yielding the final expression, 



1 ( 2 2 1 , d 2 d\ d d 



W v (x 1 ,p 1 ,x 2 ,P2)-(6) 



This is a lengthy expression, but the appearance of terms 
up to only second order together with the Gaussian func- 
tion ensures that the result, as in the case of a one-photon 
detection, is of the form of a Gaussian function multiplied 
with a second order polynomial with no first order terms 
in X2,P2- In practice, the expressions become unwieldy, 
but the analytical toolboxes in Mathematica, Matlab, or 
similar solvers readily deal with the expressions. 



IV. RESULTS FOR A REALISTIC 
EXPERIMENT 



In this section, we consider an optical parametric os- 
cillator (OPO) with mirrors with leakage rates 71, 72 and 
a nonlinear gain e. The output field from this device has 
the correlation functions 

-vi A 2 — ii 2 p -\\t-t'\ 



(a f (t)a(t') = 



71+72 4 v 2fi 2A 
71 A 2 - (i 2 r e~^ t ~ t '^ e - A l t -*' 



71+72 4 ^ 2(i 



2A 



) (7) 



where A = 2"(7i+72)+£- M = 5(71+72)— £, and where the 
time dependent field operators are normalized according 
to delta-correlated commutators [a(t), at(t')] = S(t — t'). 

The trigger mode 1 and output mode 2 are single 
modes of the field, which are in the most general case 
written as in Eq.Q with contributions from the cw 
source and from vacuum noise sources. Let us comment 
on the explicit form of these expressions: In our example 
we imagine that part of the OPO beam is transmitted 
with amplitude t through a beam splitter and the time 
dependent annihilation operator for the transmitted field 
can be written 



a T (t) = ra(t) + pb vac (t), 



(8) 



where pb vac (t) is a vacuum noise contribution, needed for 
unitarity, but causing no complications in the following. 
Field transmission losses and a finite detector efficiency 
for the trigger mode can be modelled by the effect of 
yet another beam splitter that transmits only a certain 
fraction of the field to a perfect counter and adds extra 
vacuum contributions, i.e., by a simple reduction of the 
factor t in Eq.© and by replacement of p by a corre- 
spondingly larger number. 

Further manipulation of the field prior to detection 
may be relevant, for example a frequency filter may be 



installed to avoid contributions from other modes of the 
OPO, not accounted for by Eq.Q, see for example |2lj . 
A simple frequency filter of width 7 has an exponential 
temporal transmission function, and the time dependent 
annihilation operator for the filtered field can be written 



a T {t)= T1 / expi-jit-t'^aWdt' + fabi^acit), (9) 
J — 00 

where f3ibx tVac (t) is again a short hand for the vacuum 
noise contribution. 

Eq.@ still represents a cw field, but we now assume 
that this field is being observed by some photon count- 
ing device, and that we can define sharply peaked mode 
functions centered at different times with corresponding 
annihilation and creation operators with standard single 
mode commutators. Three different counting operations 
were described in the previous section, and in a typical 
experimental implementation of the click detection the 
field will be monitored continuously in time, and even- 
tually, after an interval of no clicks, the counting device 
reports a click at time t c . We shall hence focus our atten- 
tion on the correlations between the particular discrete 
trigger mode where this happens and the remaining out- 
put from the light source. We thus define the single mode 
operator a± = J ht c (t)a T (t)dt, where ht c (t) is a normal- 
ized sharply peaked function around t = t c . If ht c (t) 
is much narrower than the temporal correlations of the 
source field, its actual shape is insignificant, and in a nu- 
merical analysis, one may discretize time in intervals of 
length dt and simply take h tc .(t) = 1/ydt on the time 
interval including t c . We note that d\ is now fully deter- 
mined as a specific integral over the source output field 
annihilation operator and vacuum contributions, defining 
thus precisely the function in Eq.QJ. 

The whole purpose of the physical setup is to pro- 
duce a useful non-Gausssian output state, and the out- 
put field, which is the component reflected by the first 
beam splitter in Fig. I, should now be made subject to 
investigation, conditioned on the click event on the trig- 
ger mode. It should be expected, that our non-classical 
output state predominantly occupies a temporal mode 
of the field that left the OPO cavity at the same time 
as the field component that caused the detector click 
on the trigger mode, and we will also expect that the 
mode has a temporal width related to the bandwidth 
of the OPO cavity as represented by the parameters A 
and p in Eq. 0) . Using the correlation functions Q , one 
may indeed determine the conditioned coherence function 
(a* (tc)S' (t)a,(t')a(t c )} , which in the spirit of standard 
photo detection theory describes the coherence proper- 
ties of the field, conditioned on a click event at time t c . 
This function is readily calculated for the specific corre- 
lation functions of the OPO, and for a low flux (small e) 
one indeed finds that this function factors into a product 
U*(t)u(t') with a mode function u(t) falling off symmet- 
rically and exponentially as function of |i — t e \, as illus- 
trated in Fig. 2a. For larger fluxes, the function no longer 
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factors, pointing to the multi mode character of the field, 
which we can understand as contributions to the field 
which are correlated with the not so distant click events 
at other times, as illustrated in Fig. 2b. 

Our formalism permits the definition of mode functions 
u(t) of arbitrary shape, and in a homodyne detection set- 
up, the extraction of a the field quadrature distribution 
of a particular mode is implemented by the correspond- 
ing temporal modulation of the local oscillator field. We 
can model losses and imperfect detectors by introducing 
the effects of an equivalent beam splitter, and hence the 
single output field mode 2 is also described by Eq.QJ 
with an explicit expression given for fatf). 

With explicit expressions for the functions fi(t') and 
the correlations (J7J|, we determine by straightforward in- 
tegration the correlations between the two discrete modes 
1 and 2, and subsequently their covariance matrix V. 
Losses and possible incoherent stray field components 
which may lead to false counts can also be included ef- 
fectively in the formalism by a modification of the co- 
variance matrix V — > LVL + N, where a diagonal matrix 
L = diag(^l - 771, VI - Vi, V 1 - W2, VI — V?) describes 
the transmission (or detection) efficiencies 1 — 77^ for the 
modes, and where N — diag(?7i+£i, r?i+£i, 772+62, 772+^2) 
describes the corresponding vacuum (rji) and additional 
(&) noise contributions, see for example Ref . \Vj\ . 

We shall now show a few examples of Wigner functions 
for output fields, obtained for different parameters. In 
the following we show the result for an OPO with one 
perfectly reflecting mirror 72 = and we shall provide 
the coupling strength of the squeezing Hamiltonian e, 
and the rates and temporal variation, characterizing the 
mode functions fi(t), in unit of 71. 

We take first a weakly excited OPO with e = O.OI71. 
We assume the triggering beamsplitter to have a trans- 
mission of 1 %, and we frequency filter the trigger field 
with a width of 7 = 571 (has only little effect on the 
field) before the detection in a narrow (jidt = 0.02) time 
window. For simplicity we assume an output mode func- 
tion of the form u(t) oc exp(— a\t — t c \, with a = ji/2. 
No losses are assumed. The Wignerfunction is shown 
in Fig. 3. It looks very much like the Wigner function 
of a one-photon number state, and we indeed compute 
the overlap with the one photon state to be as high as 
98.82 %. In the lower panel of Fig. 3. we show the results 
obtained if the output light mode suffers a 25 % loss. 
The overlap with the one photon state is now reduced to 
74.14 %. The Wigner function is negative at the origin, 
with a value of -0.3116 (no losses) and -0.154 (25 % loss), 
compared to the ideal value — 1/tt = —0.3183. 

At higher power levels, the OPO does not produce a 
one-photon state in coincidence with the click event on 
the trigger mode, but a higher excited non-classical state. 
In Fig. 4, we show Wigner functions obtained with the 
same parameters as in Fig. 3, but with e = O.271. These 
figures do not look like one-photon state Wigner func- 
tions, but they would be very attractive to produce and 
study for their non-classical behavior, and in particu- 



0.2 







-0.2 



-0.4 



0.2 







-0.2 



-0.4 




Figure 3: Upper panel: Wigner function for the single mode 
output field with mode function exp(— 71 |t— t c \/2) conditioned 
on a detector click in the trigger mode. The field comes from 
an OPO with e = O.OI71, 72 = 0, and only 1 % of the light 
is transmitted to the trigger detector. The lower panel shows 
the Wigner function obtained with the same parameters as in 
the upper panel but with 25 % loss in the output field. 



lar for their double-peaked Schrodinger cat-like character 
and for the negative values attained at the origin. With- 
out losses, we find the value of the Wigner function at 
the origin to be —0.2499 in the upper panel and —0.0889 
in the lower panel. 

We have performed the calculations with the same pa- 
rameters as in Figs. 4, but allowing for a simple variation 
of the output mode function applied in the analysis. If 
we choose a modified value of a — 0.337 instead of 0.5 for 
the decay rate of the exponential output mode function 
exp(— a\t — t c \), the negativity at the origin gets slightly 
deeper (-0.2618) than in the upper panel of Fig. 4. There 
is a host of parameters that can be varied more or less 
systematically, and if one wishes to produce a specific 
non-classical state with high fidelity, the present analysis 
is a good starting point. The natural next step is then 
to take the mode functions fi(t) as variational functions, 
and to optimize the quantity of interest (one-photon over- 
lap, negativity at the origin, cat-like behavior, ...) with 
respect to these functions. 
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Figure 4: Upper panel: Wigner function for the single mode 
output field with mode function exp(— 71 \t— t c \/2) conditioned 
on a detector click in the trigger mode. The field comes from 
an OPO with e = O.271, 72 = 0, and only f % of the light 
is transmitted to the trigger detector. The lower panel shows 
the Wigner function obtained with the same parameters as in 
the upper panel but with 25 % loss in the output field. 



V. DISCUSSION 

In conclusion we have established a formalism and we 
have presented results for the output of cw Gaussian 
light sources subject to non-Gaussian discrete measure- 
ment events on part of the field. The starting point of 
the analysis is the covariance matrix of the full Gaus- 
sian multi-mode field, given explicitly by the continuous 
time correlation functions J7J). The relevant two modes 
are described by a reduced quantum state, where one 
should formally trace over all other unobserved modes. 
The infinite dimensional covariance matrix is then re- 
duced to a covariance matrix for the two modes of in- 
terest, and rather than taking the trace, we identify this 
matrix by calculating explicitly the covariance matrix el- 
ements. The Gaussian character of the two-mode state 
is broken by the measurement, but at this stage, the sys- 
tem is so simple that any realistic trigger mode detec- 
tion scheme is readily implemented as a formal operation 
on the two-mode Wigner function of the joint system. 



This is illustrated by the action of three different non- 
Gaussian counting operations, which can all be handled 
analytically due to the simplicity of Gaussian integrals of 
polynomials. A numerical two-dimensional integral can 
of course be applied, if need be. 

We note that if an experiment is continuously not 
showing any counting events until time t c , these null 
measurements also have implications on the output state, 
and one might question the use of the steady state two- 
time correlation functions in the determination of the 
output field. In order to rely on steady state correlation 
functions, one precisely has to run the experiment con- 
tinuously in time and note that click events will occur 
after sequences of no-click intervals of varying duration, 
and that averaging over all click events corresponds to 
an averaging over all such possible histories, which, in 
turn, is equivalent to the steady state average of the un- 
probed system. This justifies the theoretical treatment, 
but it also suggests, that one may carry out experiments, 
where the output field is conditioned on a more detailed 
requirement of the trigger history, for example, that no 
clicks have occurred for a certain time before t c . This sce- 
nario will be of relevance if the trigger mode has a high 
photon detection efficiency (otherwise we automatically 
ignore many photons and hence de facto work under av- 
erage steady state conditions) . Since the no-click event is 
a projection on the vacuum state and is hence a Gaussian 
operation, we may in fact provide a simple theory also 
for this scenario. 

A useful feature of the Gaussian states is the ease at 
which the number of modes can be increased without 
rendering the theoretical treatment untractable. This im- 
plies that the above analysis readily allows generalization 
to situations where fields from several Gaussian cw and 
pulsed sources are included. We imagine in particular, 
that a recent proposal |l4| for production of one-photon 
states with no pollution from higher number state compo- 
nents involving two relatively intense OPO outputs and a 
coincidence of two single photon detection events, can be 
studied under true cw conditions. Such an analysis would 
imply a second step away from Gaussian states, and it 
remains to be seen which tools will be optimum for the 
efficient handling of this progressively more complicated 
problem. 

Despite the fact that the quantum theory of electro- 
magnetic radiation has been well established for a very 
long time and that all aspects of light generation and de- 
tection are quantum optics text book material, a "gap" 
remains between the description of light as a field and as 
composed of photons. This gap is not due to the wave- 
particle duality in quantum theory, but is of a much more 
pragmatic origin, namely that the expansion of fields in 
terms of number states is sensible and useful only for 
a single or few modes, whereas cw fields with infinitely 
many modes are more conveniently described in terms of 
the properties of field operators, labelled by continuous 
mode arguments (typically time or frequency). 

The formal equivalence between the Schrodinger pic- 
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ture expansion on number states and the Heisenberg 
picture representation of operators is useful for calcula- 
tions in the few-mode but not in the cw case. Both the 
Schrodinger and Heisenberg formalisms can be applied 
to predict mean values and variances, but the quantum 
state reduction and dynamics following a detection event 
at time t c , is not well accounted for in the Heisenberg 
picture, and the dimensionality of the Hilbert space for 
infinitely many modes makes the Schrodinger description 
unmanageable. We have previously shown that for anal- 
yses restricted to Gaussian states and operations the full 



multi mode state reduction dynamics is indeed tractable 
[l7L 12^ . I23I ]. It is this latter approach that has been 
extended in the present work, where we explicitly re- 
duce the Gaussian field description of the full multi-mode 
problem to a Gaussian field description of just two modes, 
from which the Schrodinger representation of the state re- 
duction mechanism is readily tractable, also in the non- 
Gaussian case. 
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